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Abstract. - A well known argument in cosmology gives that the power spectrum (or structure 
function) P(k) of mass density fluctuations produced from a uniform initial state by physics 
which is causal (i.e. moves matter and momentum only up to a finite scale) has the behaviour 
P(k) oc k 4 at small k. Noting the assumption of analyticity at k — of P(k) in the standard 
derivation of this result, we introduce a class of solvable one dimensional models which allows 
us to study the relation between the behaviour of P(k) at small k and the properties of the 
probability distribution f(l) for the spatial extent I of mass and momentum conserving fluc- 
tuations. We find that the k 4 behaviour is obtained in the case that the first six moments of 
f(l) are finite. Interestingly the condition that the fluctuations be localised - taken to corre- 
spond to the convergence of the first two moments of f(l) - imposes only the weaker constraint 
P(k) oc k n with n anywhere in the range < n < 4. We interpret this result to suggest that 
the causality bound will be loosened in this way if quantum fluctuations are permitted. 



In cosmology "causality bounds" are very important in various contexts. They are limits 
which can be inferred simply from the fact, which is an intrinsic feature of Big Bang models, 
that there is a finite horizon for causal processes i.e. light can travel only a finite distance 
in the time since the Big Bang. One important example is an argument, due originally to 
Zeldovich [1], which gives a strong constraint on the power spectrum (i.e. what is usually 
called the structure function in statistical physics) describing mass fluctuations. It states that 
if fluctuations are built, starting from a uniform distribution of matter, by causal physics (i.e. 
physical processes moving matter and momentum coherently up to a maximal scale), then the 
small k form of the power spectrum is P(k) oc k A . The importance of this argument is in its 
corollaries: it implies that the spectrum considered to correctly describe the perturbations at 
very large scales observed [2] in the microwave background, P{k) ~ k, cannot be produced by 
causal physics acting prior to the time when radiation decouples from matter. And it is one 
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of the motivations for and successes of the popular "inflation" model that it can produce such 
fluctuations (by modifying the causal structure of the Big Bang model at early times) . In this 
letter we describe a study of this bound - for which there is no rigourous demonstration - from 
a purely statistical physics perspective: causality is simply a bound on the distance over which 
matter can be moved coherently in building a fluctuating mass distribution from an initially 
uniform background. We introduce a set of exactly solvable one dimensional models which 
allow us to study how the fc result depends on the precise constraint which is assumed to be 
imposed on real space fluctuations. Our central result is that the condition of localisation of 
the fluctuations (in the sense usual in statistical physics i.e. finite first and second moment) 
leads only to the weaker result P(k) oc k n with < n < 4, with the fc 4 result requiring 
additional constraints on higher moments. In the context of cosmology this suggests that the 
fc 4 bound should not be taken as valid when quantum fluctuations are taken into account. 
Let us consider first a simple derivation of this fc bound on fluctuations ( 1 ). The power 

spectrum P(k) is defined as P(k) = limy_ (00 (( |<5 p (fc)| 2 ) /V) where 



6 p (k) = I d 3 xe-* k - s 5 p (x). 



v 

and S p (x) = (p(x) — po)/po is the density fluctuation field around the mean density po. 
Assuming statistical homogeneity P(k) is also given by the Fourier transform of the two point 
correlation function £(a?) = (S p (xq)5 p (xq +x)). Expanding P{k) in powers of fc one obtains, 
assuming statistical isotropy, 

P{k) = P{k) = 

J ti(x)d 3 x - f / x 2 ^(x)d 3 x + £ / x 4 £(x)d 3 x + 0(fc 6 ) 1 1 

where the integrals are now over all space, fc = |fc| and x = \x\. The fc result is obtained if 
the zero and second moments of £(x) vanish. Substituting a Taylor expansion of S p (x) ( 2 ) to 
obtain P(k), it is easy to show that these two conditions are equivalent to the requirement 

m = v „ ra »o> = „ (2) 

,,, s ( \\r($— y) 2 S (x)S (y)d 3 xd 3 y) , v 

P"(0) oc lim ^ y> M ; pyy> — =0. (3) 

That the first of these conditions Eq.|J2J follows as a consequence of mass conservation and 
the causality constraint can be understood simply. For it is just equivalent to the condition 
that the fluctuations are sub-poissonian i.e. with a squared mass variance which grows less 
rapidly than the volume. Poisson fluctuations are an upper limit because in order to obtain 
them from an initially uniform distribution, one clearly needs to move mass randomly over an 
average scale proportional to the size of the system. The causality constraint however places 
an upper bound on this scale. 

( 1 ) We note that condensed matter also provides examples of systems where the structure factor shows this 
small fc-behaviour. In particular when a system is quenched from a homogeneous (high temperature) phase into 
a broken-symmetry (low temperature) phase, the kinetics of the system is characterized by order growth [3], 
and the structure factor associated to the field describing the order has at late stages of this process a scaling 
expression S(k,t) ~ k^ l 3 F(k/k rn (t)) where k m corresponds to the maximum of 5 at a given time t. For 
systems whose phase ordering is characterized by a scalar field which is conserved during the process, the 
scaling function F(x) goes as x 4 [4,5]. 

( 2 ) In cosmology the standard derivation of the bound (see e.g. [6-9]) uses directly such a Taylor expansion 
of S p (x), which is problematic as the coefficients in this expansion are formally divergent. 
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Fig. 1 - A schematic representation of our one dimensional models. Subsequent uncorrelated steps 
of length / with probability distribution f(l) define a fragmentation of the line. The mass in each 
segment is then aggregated to a point in its centre. Mass and centre of mass are thus explicitly 
conserved locally. 

The link between the second condition Eq. (0 and momentum conservation is less evi- 
dent. Indeed the integral is not just proportional to the variance in the momentum given by 
(J v x.yS p (x)Sp(y)d 3 xd 3 y), which, in fact, diverges faster than the integration volume V. The 
necessary finiteness of the coefficient of k 2 (if it is finite) is ensured by the cancelling pieces 
which come from the coefficients of the constant and k 2 term in the expansion of S p (k) i.e. 
from the zero and second moments of the fluctuation field [5]. 

Putting aside this remaining subtlety in the demonstration of the k 4 result, we concentrate 
here on another aspect of it: the implicit assumption which has been made that the Taylor 
expansion of P(k) at k — is well defined (at least to order fc 4 ). What is the content 
of this assumption, and is it justified here? Fourier transform theory tells us that if £(x) 
is rapidly decreasing at large x, then its Fourier transform P{k) is analytic at k = 0. In 
cosmology causality might be - and often is - taken to require that physical quantities are 
strictly uncorrelated beyond the horizon scale i.e. that any correlation function is identically 
zero beyond this scale. In this case the assumption of analyticity is thus justified. Such a 
constraint on the correlation function is, however, wholly classical. In quantum field theory 
causality is imposed as a constraint on the commutators of operators at space-like separations 
[13], and does not require that correlation functions vanish. Our approach in this letter - 
through the study of the models we will now describe - is to define the constraint from causality 
as one on the form of the permitted fluctuations in real space, and then to determine what 
this constraint implies about the behaviour of P(k) as k — > 0. 

So let us consider now the following one dimensional toy model for the generation of 
fluctuations (see FigOJ. We start with a continuous exactly uniform mass density p$. Choosing 
randomly a first point, we take successive uncorrelated steps to the right of length I, where I 
is a positive number chosen randomly at each step according to the probability density f(l). 
Doing the same on the left, we define a division of the real line into segments. We then gather 
up the mass in each segment into a point, placing it in the centre of the segment. This defines 
a mass and momentum conserving discretisation of the initial mass density. The properties 
with respect to causality in this model depend wholly on those of the probability function 
/(/), as 1/2 is the maximum distance through which mass and momentum must be moved to 
create the fluctuations. We will return after our analysis of the model to discuss the subtleties 
of this constraint. Note that the discreteness of the distribution is not important here. By 
applying a smearing (e.g. gaussian) window function one could map the distribution onto a 
continuous one. Such a smearing affects only the large k properties of P(k), and thus changes 
nothing for what concerns the results below. 

We now derive exact expressions for the power spectrum in this model, and then study the 
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small k behaviour which interests us. The mass density in the segment [0, L) of finite length 
L is 



JV 



p(x) = pa ^2 l m S(x - x m ) 



(4) 



where Xi = X)}=i h + U/% an d the U are the lengths of the N consecutive segments defining 
points in [0, L). The finiteness of the first and second moments of f(l) allows to us to apply 
the central limit theorem to infer that the number of points N in such a segment is gaussian 
distributed in the large L limit about the mean N = L/(l) with fluctuations SN ~ \fh. The 
mean density is then defined i.e. we have (p(x)) — pq. From its definition it is easy to show 
that the power spectrum for k ^ is given as P(k) = hin^^oo Pl(A;) where 



P L {k) 



1 

I 



JV 



-ikx-n 



(5) 



which can be written as 




(0) 



Since we are interested in the limit L — > oo, we can evaluate the ensemble average with the 
number of points equal to the mean value N. The first ensemble average gives N(l 2 ), while 
the second, by using the definition of Xi, can be written as 



3=2 i=l 



(7) 



where (g(l)) — J °° dlg(l)f(l) for any function g(l). Performing the sum we find the exact 
expression for the power spectrum in the limit L — > +oo 



m = H> 



2Re 



(f'(l)) 2 
1 - f(k) 



(8) 



where f(k) — J °° e lkl f(l)dl is the characteristic function for the probability distribution, 
and f'(k) its first derivative. 

We can now discriminate various different cases. First let us suppose that f(k) is an 
analytic function at k = 0. All the moments of the probability distribution then exist and one 
can write the Taylor expansion 



/(*) = E 



k j (P 



3=0 



Substituted in Eq. JSJ this gives, to leading order at small k: 

, s fc 4 f(l 2 )(l 3 ) 2 (I 6 ) (l 3 )(l 4 )\ 



(9) 



(10) 
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reproducing the standard fc 4 result, with a P(k) which is analytic at k — 0. 

The derivation of Eq. © required, however, only the convergence of the first two moments. 
Let us suppose that /(/) is such that only its first n (> 2) moments are finite i.e. that 
f(l) ~ 1" /l a+1 for large /, with a > 2 (and therefore n the integer part of a). In this case the 
small k behaviour of f(k) is given by 

/(*) = E { -^r^(p) + AkX + o(k n+1 ) (ii) 

3=0 ] ' 

for non- integer a, where 

For integer a there are logarithmic corrections to these formulae. Substituting Eq. in 
Eq. (JSJ) we find the leading small k behaviour 

P(k) = C a ^k a - 2 + 0(k a - 1 ) (13) 

for 2 < n < 6, where C a is a (positive) dimensionless constant proportional to the real part of 
A. For oo > n > 6, when at least the first six moments are defined, Eq. (|10f) again is obtained, 
but with higher order corrections which are non-analytic. The non-analyticity of P(k) at k = 
implies that the correlation function £(x), which in the case that all moments of f(l) are finite 
is rapidly decaying at large separations, shows instead an algebraic tail (£(x) oc l/x a ~ r ). 

While it is the asymptotic behaviour at small k which we concentrate on here, the expres- 
sion Eq. (JSJ gives an exact expression for the power spectrum for all k. In Fig. [5] is shown 
both the result of such an exact calculation for the case f(l) = (1 + l 2 )~ 2 (i.e. a = 3), and 
a numerical simulation of the same model. The expected small k behaviour P(k) oc k can be 
clearly observed. 

Let us consider the interpretation of these results. We have been able to derive the small 
k behaviour of the power spectrum when the first two moments of f(l) are taken to be finite. 
This requirement of the fluctuations corresponds simply to the standard definition of localised 
fluctuations [10]. They are short-range correlations with a well defined characteristic scale, 
in contrast to long-range correlations (and so called "broad" distributions [10]) which lack 
such a scale. Thus we have imposed on the fluctuations that they are effectively cut-off at 
a finite length scale, which naturally then explains the bound obtained: as noted earlier the 
condition P(0) — corresponds exactly to the requirement that fluctuations in mass in real 
space should be sub-Poissonian. When we move mass coherently only up to a finite scale, 
we expect Poissonian fluctuations to form an upper bound on what one can obtain. Note 
that this is also a result which we expect to be true independently of the spatial dimension, 
as the relation between the exponents of the power spectrum and the classification (see [11] 
for a discussion) of mass fluctuations in poissonian, sub-poissonian and super-poissonian is 
independent of dimension ( 3 ). Note also that the correlation function, as we have seen, decays 
at least as fast as a power-law with exponent greater than the dimension (a > 2), which indeed 
corresponds to an effectively short-range correlation with an integrable correlation function 
(whose integral is zero). 

( 3 )A generalization of this model to arbitary dimensions will be presented elsewhere [12]. 



(12) 
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Fig. 2 - The power spectrum measured in a simulation of our one dimensional model with f(l) — 
(1 + l 2 )" 2 (i.e. a = 3) for a box size L = 400, and one thousand realisations. The black line is 
the exact result obtained directly from Eq. JHJ. Both curves show the predicted small k behaviour 
P{k) oc k. 



While we have this simple interpretation of our results, what we have found is actually very 
non-trivial and surprising in another respect: the exponent at small k, which characterises 
the large scale fluctuations is not universal for this class of localised distributions with finite 
mean and variance. Rather the exponent depends in a subtle way on the convergence of higher 
moments. We stress that this means the small-k behaviour is thus determined by the large 
masses (i.e. the tails) of the these distributions despite the fact that they are not broad. 

Now let us return to the interpretation of these results in relation to the cosmological 
causality bound, which was our original motivation. The problem of the causality bound can 
be addressed with this simple model, because 1/2 is the maximum distance through which mass 
and momentum must be moved to create the fluctuations. If we impose that f(l) is strictly zero 
above some finite scale l c , no mass can be moved over a scale larger than this (l c is then the "size 
of the causal horizon" ) . Such a constraint is what one would take to be imposed classically 
by causality. And our model gives in this case the fc 4 result, in accord with the standard 
cosmological bound. However our model shows that a generalisation of this result to the case 
where the fluctuations are cut off at a characteristic scale l c , but smoothly, is non-trivial. And 
it is precisely such a generalisation which is needed if the k 4 bound is to extended to the case 
in which quantum fluctuations are admitted (i.e. that the physics generating the fluctuations 
is quantum in nature). This is so because quantum correlation functions generically do not 
vanish at any separation, as the required sharp localisation of states is impossible (cf. Reeh- 
Schlieder theorem [14]). This means that to extend the standard result also to include quantum 
fluctuations, one must assume that any (causal) quantum mechanical evolution giving rise to 
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fluctuations in an expanding FRW Universe (with sub-luminal, non-inflationary expansion) 
gives rise to, for example, a strictly rapid decay of correlation functions at super-horizon 
scales, and excludes the rapid power-law decay admitted by a generic definition of localisation 
of fluctuations. While a rapid decay as an exponential is indeed typical of purely quantum 
fluctuations in the simplest cases one can envisage (e.g. the quantum correlations at space- 
like separations of a free massive scalar field in its ground state) such an assumption is too 
restrictive a form of the condition that quantum fluctuations be localised. An interesting 
analogy is given by screening of electrostatic forces in a neutral plasma, which is analagous 
to the effect of the existence of a causal horizon relative to a theory in flat space, as it leads 
to a physical cut-off scale beyond which the effect of a source is no longer felt. While such 
screening gives rise to exponential decay of the (charge-charge) correlation function in the 
case of a Coulomb potential, it can be shown easily [15] that, for a 1/r 2 potential, the same 
physical mechanism leads to an effective short-range potential which decays algebraically (as 
— 1/r 4 in three dimensions). Indeed the spectrum of fluctuations (at thermal equilibrium) at 
small k in this case behaves as P(k) oc k, instead of P(k) oc k 2 in the case of electrostatics [15]. 
This suggests that a power spectrum of fluctuations P(k) oc k at small k might, for example, 
be obtained from an appropriate interacting field quantized in a (sub-luminally expanding) 
FRW background. Such a model would provide an interesting alternative to the standard 
(non-causal) inflationary mechanism of producing such fluctuations at superhorizon scales. 
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